The Giant Dipole Resonance as a quantitative constraint on the symmetry energy 
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The possible constraints on the poorly determined symmetry part of the effective nuclear Hamil- 
tonians or effective energy functionals, i.e., the so-called symmetry energy S{p), axe very much under 
debate. In the present work, we show that the value of the symmetry energy associated with Skyrme 
functionals, at densities p around 0.1 fm~'^, is strongly correlated with the value of the centroid of 
the Giant Dipole Resonance (GDR) in spherical nuclei. Consequently, the experimental value of 
the GDR in, e.g., ^"^Pb can be used as a constraint on the symmetry energy, leading to 23.3 MeV 
< S{p = 0.1 fm-^') < 24.9 MeV. 



oo 

o 

O 

(N 

X) 
CD 

vn 



> 

00 

\o 
cn 

(N 
O 
OO 

O 



I. INTRODUCTION 

The nuclear structure community is currently striving 
to determine a nuclear energy functional as universal and 
as accurate as possibile. Extraction of this functional 
from a more fundamental theory like QCD is of course 
desirable, and there has been progress along this line. 
At present, however, it is still unavoidable to work with 
functionals which depend on free parameters which have 
to be determined by some fitting procedure. 

Existing functionals include those based on a covariant 
formulation [1] as well as those based on nonrelativistic 
formulations. Restricting ourselves to the latter case, we 
note that a nuclear energy functional can be defined in 
a general way, without being derived from an underlying 
Hamiltonian in conjunction with a reference state. How- 
ever, most of the existing functionals to date are derived 
from an effective Hamiltonian Hcs which includes the ki- 
netic energy plus a two-body interaction. In this case, the 
total energy is the expectation value of Hcs over the most 
general Slater determinant |(f>). Both zero-range interac- 
tions like the one proposed by Skyrme at the end of the 
fifties, and systematically parametrized for Hartree-Fock 
(HF) calculations since the seventies or finite-range 

interactions like the Gogny force [3| , lead to satisfactory 
descriptions of many nuclear properties. 

Our tool of choice in the present work is the zero-range 
Skyrme force, from which one can derive a functional £[p\ 
which is a function of local densities only. For a system 
that is not symmetric in neutrons and protons, the total 
energy depends both on neutron and proton density: 



E[p] 



ds.r £ip„ir),pp{r)). 



(1) 



For the sake of simplicity, we have not indicated that 
in general the energy depends not only on the spatial 
densities, but also on gradients Vpg, on the kinetic energy 
densities Tq and on the spin-orbit densities Jq (where q 
labels n,p) [1, @. 

In infinite matter, one has a simple expression in terms 
of the spatial densities only. Instead of p„ and pp, one 
can use the total density p and the local neutron-proton 
asymmetry. 



This quantity should not be confused with the global 
asymmetry (N — Z)/A. In asymmetric matter, we can 
make a further simplification on f (p, S) by making a Tay- 
lor expansion in S and retaining only the quadratic term, 



Sip,S) 



£q{p,S 
£o{p,S 



0)+£sym(p)(52 

0) + pSip)5'. 



(3) 



The first term on the r.h.s. is the energy density of sym- 
metric nuclear matter fnm , while the second term defines 
the main object of the present study, namely the sym- 
metry energy S{p). The symmetry energy at saturation 
S{pq) is denoted by different symbols in the literature: 
J, Or or 04. We stress that Eq. ([3]) is not really a sim- 
plification: the coefficient of the term in S"^ which should 
follow, for the Skyrme parameter sets employed in this 
work, is negligible at densities of the order of po. We 
remind that the pressure of the system can be written in 
a uniform system as 



P 



dE 
dV 



dp p 



(4) 



Pn - Pp 



(2) 



This quantity is evidently related to the density depen- 
dence of the energy functional and of the associated sym- 
metry part defined above. 

The magnitude and the density dependence of the sym- 
metry energy S{p) are not yet well understood 0|- In 
brief, there exist at present three main research lines 
aimed at constraining the behavior of the symmetry en- 
ergy, by using either nuclear structure data, or observ- 
ables related to heavy-ion collisions, or evidences from 
the study of neutron stars. 

Within the realm of nuclear structure, the symmetry 
energy affects of course all properties of nuclei having 
neutron excess, including basic ones like masses and radii. 
In particular, much attention has been focused on radii 
since Typel and Brown [8j, |9|] have noted that the neu- 
tron skin thickness SR = (r?V/^ — {rp)^^'^ is correlated 
with Pnm(p = 0.1) (see also The issue is also inves- 

tigated in Ref. 11], where the correlations between the 
neutron skin thickness and other quantities are discussed. 
The experimental accuracy is not sufficient (so far) to 
limit the acceptable range of the neutron skin thickness 
so that this can constrain a given equation of state; the 
Parity Radius Experiment (PREX) at JLAB promises to 
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achieve this task [T3| ■ Another interesting way to fix the 
value of the neutron skin thickness, and extract informa- 
tion on the symmetry energy, is to go through the study 
of the isovector spin-dipole resonance (SDR) sum rule; 
also in this case, the experimental difficulties hinder a 
too definite conclusion (see Ref. [13]). 

In this paper we shall instead concentrate on the cor- 
relation between the symmetry energy and the energy of 
the Giant Dipole Resonance (GDR). The idea is not new, 
but in the present work we develop it based on a fully mi- 
croscopic approach, namely within a self-consistent Ran- 
dom Phase Approximation (RPA) scheme to calculate 
the GDR properties. In the past [13], as well as in re- 
cent works [151 , the connection with the symmetry energy 
has been discussed starting from a macroscopic, hydro- 
dynamical description of the GDR (in particular, using 
the Steinwedel- Jensen ansatz which is known to be not 
fully reliable). Consequently, we believe that our results 
are more relevant from a quantitative point of view. 

If one tries to constrain the symmetry energy by means 
of the study of heavy-ion collisions, or by neutron star 
observables, it is likely that somewhat different physics 
is involved. In general, the behavior of the symmetry 
energy on a broader range of densities is involved. In 
heavy-ion collisions maximum densities up to ^4-5 pa 
can be attained [l^. On the other hand, data at lower 
incident energies are believed to be able to constrain the 
nuclear EOS below po [13] • The study of neutron stars, as 
well, brings in the physics of both low-density and high- 
density neutron matter (for a comprehensive review, cf. 
Ref. [18]). A recent study [l^ which is close in spirit to 
ours, has examined a large set of Skyrme forces, trying 
to determine those which have a satisfactory behavior 
in reproducing the neutron-star observables. One should 
remark, however, that there are many caveats in the lit- 
erature against the use of functionals in a density regime 
far from that in which the functionals are fitted and usu- 
ally employed. In particular, we mention here that Monte 
Carlo calculations of neutron matter at low density (20| 
show that in this regime E/A is about one half of the 
Fermi energy of the non- interacting neutron gas, and this 
behavior is not reproduced by any effective mean field 
functional. 

Consequently, we do not discuss in detail in the present 
work the possibility of an overall constraint on the sym- 
metry energy extracted by different kinds of studies. We 
briefly discuss in the conclusions to what extent our re- 
sults can be compared with a few others in the literature. 



II. THE CORRELATION BETWEEN THE GDR 
AND THE SYMMETRY ENERGY 

Our starting point will be the hydrodynamical 
model of giant resonances, proposed by Lipparini and 
Stringari [2l| . They assume an energy functional which is 
simplified yet sufficiently realistic, solve the macroscopic 
equations for the densities and currents, and extract ex- 



pressions for the moments mi and to_i associated with 
an external operator F [ruk = / dES{E)E^ where S is 
the strength function associated with F) . The expression 
for mi is proportional to (1 -|- k), where k is the well- 
known "enhancement factor" which in the case of Skyrme 
forces is associated with their velocity dependence. The 
expression for m_i, in the case of an isovector external 
operator, includes integrals involving £sym and F. They 
can be evaluated in a simple way if one assumes the valid- 
ity of the leptodermous expansion. We write the volume 
and surface coefficients of the expansion of Esym as 6voi 
and 6surf, respectively. By specializing F to the isovec- 
tor dipole case, the following expression is obtained (for 
details, cf. Ref. (U) 



mi 
m_i 



\ 



Ovol 



m{r'^) 



O Ovol 



(I + k). (5) 



This equation yields values of the centroid energy which 
are in rather good agreement with those of microscopic 
RPA calculations. It turned out to be useful in a previous 
study [2^ . in order to constrain directly the parameters 
of the isovector part of the Skyrme interaction. Here 
we shall use it as a guideline, and try instead to find a 
quantitative connection between the energy of the GDR 
and the symmetry energy. 

The ratio can be evaluated through the calcula- 
tion of a semi-infinite nuclear slab. This has been done, 
e.g., in Ref. [11] (cf. their Sec. 3.2.3). We do not discuss 
here the approximations made in the derivation, but we 
use the fact that the mentioned ratio can be written in 
terms of the symmetry energy and its derivatives. If we 
insert this result into Eq. ([5]) we obtain 



E^i = 



m{r'^) 



where 



9Aip) 



5a(po)(1 + k), 



S{p) 



(6) 



(7) 



For a given heavy nucleus we can safely consider the first 
of the three factors under the square root as a constant, 
since different Skyrme forces do not vary widely in their 
predictions for (r^). We have evaluated the term gA{p) at 
p — Po for v4=40, 124, 208, and for a number of Skyrme 
forces. We have found that it is strongly correlated with 
the value of S{p) in the range p = 0.08 - 0.12 fm"^. 
The specific case ^=208 is displayed in Fig. [T] in the 
case p = 0.1 fm~"^, for which the correlation coefficient is 
maximum. 

Although we have not been able to deduce this correla- 
tion in an analytic way from the expression of the Skyrme 
functional, this result, together with Eq. ([5]), motivates 
us to look for a direct correlation between the centroid of 
the GDR and the symmetry energy, through the quantity 



/(p) ^ ^S{p){l + n), 



(8) 
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for p ~ 0.1 fm""^. In the next Section, we discuss this 
correlation for p — 0.1 im~^ in the case of ^"^Pb. We 
also discuss in some detail the choice of the Skyrme forces 
that we have employed. 
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FIG. 1: Correlation between the quantity (7a=208(po) and the 
symmetry energy S'(p) for the Skyrme forces listed in Table U 
at p = 0.1 fm~'^. The value of the correlation coefficient is 
r = 0.981. 



III. RESULTS 

We have obtained results for the GDR in ^"*Pb by us- 
ing a series of microscopic Hartree-Fock (HF) plus Ran- 
dom Phase Approximation (RPA) calculations. Skyrme- 
RPA theory is well known since many years, especially 
in its matrix formulation. Recently, we have developed 
a scheme which is fully self-consistent and is discussed 
in [2^ . There is no approximation in the residual inter- 
action, in that all its terms are taken into account (includ- 
ing the two-body spin-orbit and Coulomb interactions). 
The occupied states are determined by solving the HF 
equations in a radial mesh extending up to 24 fm. The 
continuum is discretized since the positive energy states 
are obtained using box boundary conditions. Particle- 
hole (p-h) configurations which constitute the basis for 
the RPA matrix equations are included up to typically 
60 MeV so that the value of mi is at least 98% of the 
well-known value obtained from the double commutator. 
In a few cases, due to the instability of RPA, we had to 
resort to Tamm-Dancoff Approximation (TDA) calcula- 
tions but this does not affect significantly our results. 

The calculations have been performed for a set of 20 
Skyrme interactions. It is well known that more than 100 
Skyrme parametrizations have been proposed in the lit- 
erature, including some which have been employed only 
in limited and specific cases. It is hard to define in a 
clear-cut way a "standard" subset to be analyzed; how- 
ever, we have decided, in the present context, not to con- 
sider parametrizations which (a) have an associated K^a 



outside of the range 210-270 MeV, in keeping with the 
conclusions reached in Ref. [HI by studying the giant 
monopole resonance (GMR) in ^"^Pb, (b) reproduce the 
experimental value of the GDR in ^o^Pb (13.46 MeV [111) 
within ± 2 MeV. 

We have determined our set including forces proposed 
by different groups, and at different times. In this sense 
the set can be considered representative enough. In the 
cases in which several forces have been proposed in the 
same reference, we have not included more than two 
forces, in order to avoid a too strong bias. The forces 
are listed in Table [H we also provide the reference from 
which the parameter sets have been taken, the value for 
E^i{RPA) EE obtained from our RPA calcula- 

tion, and the values of the quantities /(O.l), S'(O.l) and 
K (cf. Eq. ®). 





Ref. 


E-i 


/(O.l) 


5(0.1) 


K 






[MeV] 


[MeV'/''] 


[MeV] 




SkA 


[23 


15.14 


6.27 


23.43 


0.68 


SkM 


M 


13.94 


5.65 


23.26 


0.37 


SGI 


128] 


14.16 


5.62 


20.40 


0.55 


SGII 


ggj 


13.56 


5.34 


20.98 


0.36 


SkM* 


[29] 


13.89 


5.64 


23.26 


0.37 


RATP 


[30] 


15.17 


6.06 


23.55 


0.56 


SkT4 




11.47 


4.86 


23.58 


0.00 


SkT6 


m 


12.17 


4.92 


24.20 


0.00 


Rs 


[32] 


12.76 


5.20 


20.05 


0.35 


Gs 


[32] 


12.62 


5.20 


20.03 


0.35 


SkI2 


[33J 


12.29 


5.00 


21.17 


0.18 


SLy230a 


m 


12.49 


5.04 


25.43 


~ 


SLy4 




13.40 


5.45 


25.15 


0.18 


SLy5 


[5] 


13.28 


5.42 


24.88 


0.18 


SkO' 


[35] 


13.85 


5.00 


22.61 


0.11 


MSk7 


[36] 


12.10 


4.86 


24.56 


-0.04 


vllO 


[37] 


12.13 


4.80 


24.79 


-0.07 


v075 


[37] 


13.97 


5.62 


25.30 


0.25 


SK255 


[38J 


13.98 


5.94 


25.39 


0.39 


LNS 


[39] 


13.95 


5.43 


23.20 


0.27 



TABLE I: For the Skyrme parameter sets considered in this 
work, we provide the values of E-i{RPA), /(O.l), 5(0.1) and 
K. All these quantities are defined in the text. 

We find a strong linear correlation between the val- 
ues of E-i{RPA) and /(O.l), which are shown in Fig. [2] 
together with the interpolating straight line /(O.l) = 
a + bE^i{RPA). The value of the correlation coefficient 
is r = 0.909. Before discussing the extraction of the value 
of the symmetry energy, we should stress that we have 
not been able to correlate the GDR simply with S{po); 
this may be possible (cf., e.g., Ref. [13]) at the price of 
restricting oneself to a small set of Skyrme forces. 

We can now make avail of the experimental values 
of the GDR centroid E^i{exp) and of the enhancement 
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factor K to deduce the best value of the symmetry en- 
ergy. While the value of E_i{exp) in 208p|^ ^^as been 
rather well determined from photoabsorbtion measure- 
ments, E-i{exp) — 13.46 MeV [l^, the value of k is less 
precise. Since 

J a{E)dE = 60 ( — j •(! + '«) MeV • mb, (9) 

K would be determined if the integrated photoabsorb- 
tion cross section <t{E) had been mesured up to large 
energies (essentially, up to the pion production thresh- 
old). For obvious experimental difficulties, the photoab- 
sorbtion cross section has been measured up to much 
lower energies (26.4 MeV in the case at hand, namely 
^"^Pb [4]|). The uncertainty on k has been already esti- 
mated in [21[, where it has been stated that k should lie, 
approximately, between 0.2 and 0.3. We have checked 
this in some detail in the following way. From the ex- 
perimental cross section measured up to the maximum 
energy of 26.4 MeV, that is, from Eq. ([9]) replacing the 
upper limit in the integral with 26.4 MeV, we have de- 
duced an "elfective" k, that we can call ^26. 4- The Skyrme 
forces which have an associated K26.4 within the experi- 
mental limit are those which have k, between Kmin^O.lS 
and «;max=0.26. In the following, we have used these val- 
ues. Our best value for k is of course an average between 

^min and /i^max- 




13 14 
E , [MeV] 



FIG. 2: Correlation between the energy of the GDR and the 
quantity /(O.l). The definition of these quantities, and the 
related discussion, can be found in the text. 



At this point, the best value for 5(0.1) is found as 

(a + bE!_Y)^ 



3(0.1) = 



(10) 



(l + «) ' 

where a, b come from the fit. The error is obtained from 



exp\2 



(11) 



where at comes from the fit and the variance is 
calculated with respect to the interpolating straight line. 
Having determined the ±1<t interval around the mean 
value for the quantity S'(0.1)(l + k), we obtain 



5.419 r—— 5.422 

< ^/S{0.1) < 



(12) 



This can be considered one of the main results of the 
present investigation. By introducing the values of Kmin 
and Kmax discussed above, one obtains a further (more 
direct) constrain, that is. 



23.3 MeV < S'(O.l) < 24.9 MeV. 



IV. CONCLUSIONS 



(13) 



Using a representative set of Skyrme effective func- 
tionals we have found a clear correlation between the 
energy of the GDR in ^"^Pb and a simple function of 
the symmetry energy at density p ~ 0.1 fm~^ and of the 
enhancement factor k associated with the velocity depen- 
dence (and with the effective mass) of the various func- 
tional. Using the well established experimental value of 
the GDR we have extracted a range of acceptable values 
of for 5(0.1) (cf. Eq. ([T3D). 

It would be important to test whether other classes of 
effective functionals lead to a similar result. More gen- 
erally, it will be essential to study the interplay between 
constraints coming from the different kinds of works men- 
tioned in the Introduction, which deal with different en- 
ergy and density regimes. In fact, a better knowledge 
of the symmetry part of the nuclear effective function- 
als, and in particular of its density dependence, would 
be highly instrumental for the study of systems ranging 
from exotic nuclei to pure neutron matter. 
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